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Abstract 

The conventional theory of resonance broadening for a two-species plasma in a mag- 
netic field is revised, and applied to an ionospheric turbulence case. The assumptions 
made in the conventional theory of resonance broadening have, in the past, led to 
replacing the frequency u by u + ik\D* in the resonant part of the linear dielec- 
tric function to obtain the nonlinear dielectric function. Where D* is an anomalous 
diffusion coefficient due solely to wave scattering of the particle orbits. We show 
that in general these assumptions are not valid, and consequently the straightforward 
substitution of frequencies is not legitimate. We remedy these problems and derive 
expressions for the time-dependent components of the diffusion tensor. The improved 
resonance broadening theory is developed in the context of an ionospheric problem, 
namely that of the Farley-Buneman turbulence in the auroral E-region. A kinetic 
description of the electrons is used. A general expression for the nonlinear dielectric 
function is derived in the special case where no parallel electric field is present, and 
the differences with the conventional dispersion relation are discussed. 

1 Introduction 

The problem of wave-particle interaction has always played a critical role when one tries 
to understand the saturation mechanisms for plasma instabilities. Linear theory, which is 
a single wave theory, which does not take into account the wave-particle interaction fails 
to conserve energy and momentum, and consequently does not and should not predict 
saturation. This is unphysical. But it is important to remember that the linear theory 
ceases to be valid after a trapping time; the time neccessary to a charged particle to bo- 
unce near the bottom of a potential well of a finite-amplitude wave. Consequently, the 
predictions of linear theory beyond the trapping time are not valid. A first remedy to this 
fundamental problem is provided by the so called Quasilinear theory, a weak turbulence 
theory which takes into consideration the wave-particle interaction and requires a slow 
time dependence of the background particle distribution function. The ensemble averaged 
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distribution function is then a solution to a Fokker-Planck equation. In the absence of 
sources and sinks, quasilinear theory is then described through a diffusion equation which 
predicts saturation of the instabilities when the background distribution function becomes 
constant, or "plateaus", along the diffusion paths which are the characteristics of the par- 
tial differential equation. However, in the presence of sources and sinks, such as a back- 
ground electric field or collisions, there is no saturation, since the sources and sinks tend to 
destroy the plateau, and consequently the unstable waves keep on growing until they be- 
come large enough for the other nonlinear processes to enter the picture. This is the case 
of the problem we have elected to address, namely that of the modified two stream Farley- 
Buneman instability, where the source of free energy is the electrojet and the sink is due 
to the collisions of the electrons with the background neutrals. In other words quasilinear 
theory does not provide an ultimate saturation mechanism much needed to describe steady 
state turbulence and predict reasonable saturation amplitudes. Nevertheless the quasi- 
linear theory has stood the test of time when it comes to predicting the onset of instabilities 
to a certain degree of accuracy. A first attempt to improve on the quasilinear theory was 
described by Dupree [5] and Weinstock [20] for an unmagnetized plasma and by Dupree [6] 
and Dum and Dupree [7] for a magnetized plasma. Dupree [5] derived a perturbation 
theory based on the knowledge of the electric field which allowed him to find the exact 
particle orbits which were then used in the perturbation solutions of the Vlasov equation. 
The principal result of this improved perturbation theory is a broadening of the wave- 
particle resonance function, which in the conventional quasilinear theory is S(u> — k • v). 
This result has been the subject of many debates on strong turbulence theories in the 
past two decades. The validity of some of the assumptions made by Dupree [5] was 
questioned by a number of authors (Cook and Sanderson [3], Rolland [15]). A number of 
authors (Salat [18], Ishihara and Hirose [10], Ishihara et al. [11, 12] have since taken the 
task to study the problem of resonance broadening very thoroughly, and were able to show 
successfully the shortcomings of the conventional resonance broadening approach. They all 
addressed the resonance broadening problem in the case of unmagnetized plasmas, as well 
as the case related to drift waves in a shear magnetic field because of the mathematical 
tractability and its reduction to a form similar to the unmagnetized case. Kleva and 
Drake [21] addressed the problem of stochastic E x B particle transport using a dynamical 
system approach. More recently Kleva [14] investigated the problem of energy transport 
in a magnetically confined plasma by E x B flow generated by a spectrum of electrostatic 
waves stil using a dynamical system's approach. He found like Ishihara et al. [12] that 
the diffusion coefficient scales like .E 4 / 3 instead of E 2 as predicted by the conventional 
quasilinear theory. The problem addressed by Dupree [6] and Dum and Dupree [7], that 
of the broadening of the wave-particle resonance in the presence of background magnetic 
field, has not been addressed since probably because of its complexity. 

In this paper we have investigated the problem of resonance broadening for a turbulent 
magnetized plasma and compared our results to the classical calculation of Dum and 
Dupree [7]. We have investigated the implications of the improved results on the Farley- 
Buneman instability that occurs in ionospheric plasmas, and shown that one cannot obtain 
the nonlinear dielectric function by just substituting the frequency w by u> + ik\D* in the 
resonant part of the dielectric function. In Section 2 we describe the mathematical model 
used to investigate the Farley-Buneman instability, and derive the different components of 
the diffusion tensor. In Section 3 we derived the generalized dispersion relation, and finally 
in Section 4 we properly reinterpret the thresholds conditions for the Farley-Buneman 
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instability in the absence of a parallel electric field. 

2 The Mathematical Model 

2.1 The Quasi-Linear Approximation 
2.1.1 The Ion Description 

The ions are assumed to be highly collisional in the region of interest, namely the auroral 
E-region, and unmagnetized. The ion convection is also assumed to be negligible, i.e., the 
nonlinear ion terms are neglected. 

With these assumptions a linear fluid model is adopted for the ions, and is best de- 
scribed by the linearized momentum and continuity equations. Following Sudan [19] we 
write 

m i~KT = - eV< /> - — V<5n - v in mi5vi, 
ot no 

— + n V-5v l = (2.1) 

operating on the continuity equation (2.1) with ^— + Vin^j leads to 

d \ dSn ( d \ ( e T 1 

57 + v in — = -n V • — + u in )Svi = n V • { — V<f> + ^-V5n \ . (2.2) 
ot J ot \Ot J y rrii TiQirii J 



Equation (2^) allows us to express the density fluctuation in terms of the electric field. 



Taking the Fourier transform in space and time of equation ( |2.2| ) leads to 

ji _ £2_»_ + i Vin(j j \ $ nkuj = r ^L]^ ( j )Uw (2.3) 
mi J m 

This equation will eventually be used in Poisson's equation to determine the dispersion 
relation. 

2.1.2 The Electron Description 

The electrons are described by the Vlasov equation with a relaxation model for collision 
operator. The electron distribution function satisfies the following equation 

^ + V . V - — (E+- XB) .V v ) / e (x,V,i) = 

ot m e \ c / I 

(2.4) 

-Ven (/ e (x,V,t) - - f ) , 

where 

E(x,t) =E + 5E(x,t), 

(2.5) 



n(x,t) = J dvf e (x,v,t). 
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The electric field has been separated into two parts; Eo represents the electrojet back- 
ground electric field, and <5E the fluctuating field. 

This leads to the following distribution function, which is in turn separated into a weakly 
space and tine dependent average distibution function (/ e (x, v, t)), and a fluctuating part 
5f e (x, v, t), that is 

/ e (x, V, t) = (/ e (x, V, t)) + «5/ e (x, V, t). (2.6) 

We should point out that both the electric field and the distribution function decom- 
positions are exact. One can easily write equation fl2.4|) in a compact form 

^ + i£(x, v, A / e (x, v, t) = ^E(x, t) ■ V v / e (x, v, t), (2.7) 

where the operator C is given by: 

g / v \ fo f 

i£(x, v, t) = v • V I E H — x B ) • V v + v en — v en — / dw. 

m e V c / no 7 

Note that the collision operator is a linear operator, and consequently its incorporation 
into the collisionless Vlasov equation is not as complex as if we would have used a Fokker- 
Plank or any more sofisticated collision operator. 

Equation ( |2.7j ) can be solved using different techniques. In the quasi-linear approx- 
imation one considers only the wave particle interaction and neglects any other effects, 
such as wave coupling and radiation effects which are nonlinear effects of higher order in 
a perturbation analysis based on the amplitude of the wave electric field. 

To solve the quasi-linear problem we start by ensemble averaging the collisional Valsov 
equation to obtain 

d + Ven + V • V - — (Eo + - X b) • V v J (/ e ) 



Of 



(2i 

Venh + — Vv (SESf e ). 



The next step consists of substituting the expression for the distribution function into 
the collisional Vlasov equation to obtain 



^ + v en + v • V - ^ (e + \ x B)- V v ) ((f e ) + 6f e ) - ^5E ■ V v (./; ) 



m, P \ c J J m, 



fo & 

Venh + V en Sn H V v • (SE5f e 

n m e 



(2.9) 



where 8n is defined as 
5n = dv5f e . 



Taking the difference between equation (|2.9|) and equation (|2.8D leads to 

^ + v en + v • V - — (e + - X b) • VvWe = 
at m e \ c J J 

-—SB ■ V v (/ e ) - VeJ—Sn = — V v ■ {SB5f e - (5B5f e )} . 



(2.10) 
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The quasi-linear approximation consists of neglecting the right hand side of equa- 
tion ( p. 10 ) which describes the nonlinear mode coupling terms. When this term is ne- 
glected, the equation to lowest order in the electric field becomes 

<} ■ - + v • V - — (e + - x b) • V v ) 8f e = v en —6n + — 5E • V v (/e). 



ar m e \ c J J uq m, 

This equation can be solved formally using Fourier transforms and defining the following 
operator 

Gku,(v) = <^ -iuj + v en + ik ■ v - — (E + - x B ) • V v 

[ m e \ c J 

The expression for the fluctuating distribution function can then be written in the following 
form 

*Wv) = — ■ G kw V v (/ e ) + G^/o. (2.11) 

m e n 

Integrating over velocity space one obtains an expression for the density fluctuation which 
when substituted into expression ( p. 11 ) for the fluctuating distribution function leads to 

e ^ SE ^> • / dvG kw V v (/ e ) 

fifelau = <5E kw • GkcjV v (/ e ) + V en — T~J~~F\ 1 Gk^/o (2-12) 

m e n - Ven J avG^/o 

when substituting expression (|2.12j ) into the equation governing the evolution of the aver- 
age distribution equation one obtains 

d(f e ) _ d p d(f e ) d 

dt <9v = <9v dv ~ e 

In other words one obtains a Fokker-Planck equation with the diffusion and drag coeffi- 
cients defined as follows 

e 2 

D = —5 (<5E5EG) , 
— 

2 / r^vQ^iM \ 

F(/ e ) = E (/ e ) - — ^ ( 5E5EG • TTT^T ^ ) • 

m e m| \ n - ^ eri J dvG/ / 

It is clear that the drag term in the Fokker-Planck equation is solely due to the background 
electric field of the electrojet and to the electron-neutral collisions. 

In the absence of sources or collisions the conventional quasi-linear theory predicts a 
saturation of the fluctuating fields when the distribution function becomes constant along 
the diffusion paths. That is, the quasi-linear theory predicts a zero growth rate. This 
however, is valid only when there are no sources or collisions; when such effects are present, 
like in our case, the distribution function never plateaus along the diffusion paths since 
the source and the collisions tend to destroy the plateau, which in turn leads to a non zero 
growth rate; the oscillation amplitude continues to grow until nonlinear processes enter 
the picture. In other words, there is no ultimate saturation of the fluctuations through 
the quasilinear process. 



Resonance Broadening Theory of Farley-Buneman Turbulence 



443 



This leads us to conclude that the conventional quasi-linear theory is not the ultimate 
stabilization mechanism, and improvements on the theory are needed. One possible theory 
that has been suggested, as a first attempt to remedy the problem from which the quasi- 
linear theory suffers, is the resonance broadening theory which was first introduced by 
Dupree [5], and applied by a number of authors, see for example Sudan [19], Robinson [16], 
and Robinson and Honary [17], to the problem of irregularities in the E region. 

In the next section we shall develop the resonance broadening in some details. We will 
also discuss, to a certain extent, the validity limits of the conventional resonance broad- 
ening. We remedy some of the problems, and show the shortcomings of this wave particle 
interaction, as well as the need to provide for a complete theory of plasma turbulence. 
The latter will be addressed in a companion paper describing a fully selfconsistent kinetic 
theory for the Farley-Buneman instability. 



2.2 The Resonance Broadening Approximation 

In this section we shall address the problem of wave particle interaction in a magnetized 
plasma through a nonlinear formalism that includes the nonlinear effects of the waves on 
the particle orbits, but not vice versa. The electric field and the corresponding distribution 
function are decomposed according to equations (|2.5| ) and (|2,6]). 

The fundamental assumption and the goal of the resonance broadening theory is to 
evaluate the modification of the quasilinear resonance between the particles and waves. 
This theory neglects the coherent contributions which arise from the coupling between 
the waves and the background oscillations, as well as interactions between the background 
oscillations. 

We define the following operator C is then defined as follows 



i£(x,v,t) =v-V E + £E + - xB • V v + co "(x,v,f), 

m e V c J 

where the collision operator O coU is given by 

O co "(x, v,t) = v en - v e j^± [ dV 

no J 

which allows us to rewrite equation (|2.4D in a compact form 
d 



+iC(x,V,t) ] / e (x,V,t) =0. 

The equation describing the evolution of the ensemble averaged electron distribution 
function is given by 

^ + i£o(x, v, if) (/ e (x, v, t)) = -^-V v • (<5E(x, t)6f e (x, v, t)), 

where Cq is a linear operator defined by 

i£o(x, v, t) = v • V + ( — - E + fl e x v J • V v + v en - vj-^- [ dV. 

\ rn e ) n J 
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The equation governing the evolution of the fluctuating part of the electron distribution 
function can be written as follows 

J^+iC(x,v,i)) <5/ e (x,v,i) = ^-<5E(x,t) • V v </ e (x,v,t)>, (2.13) 

where as one can see the nonlinear term 6E5f e . This term will take into account the effects 
of the waves on the electrons in the resonance broadening approximation, but will not 
address or investigate the effects of the electrons on the waves, neither does it take into 
consideration the wave-wave interaction. The "Resonance Broadening" approximation 
assumes that the quadratic nonlinearity leads to a nonlinear correction to the particle 
orbits, which should be taken into account when the fields become large enough. 

Equation (|1|) can be solved using a Green's function analysis, i.e., we define the 
Green's function Q which satisfies 

^ + iC(x, v, t)\ g(x, v, t; x', v', = <5(x - x')<5(v - v')6(t - t'). (2.14) 

Following Ishihara et al. [12] the operator iC can be rewritten as follows 
iC = iCq + iCi, 

e -E + n e x v ) • V v + C coU 



iC = v-V+l Eo + flexv .v v + O co ", ( 2.15) 



id = — —5E(x,t) ■ V v . 

m e 

The electron distribution function is then obtained 

f e (x,v,t) = j dx'dVg(x,v,t;x',v',t )f e (x,v',t ). 

Before we get into the details of the derivation, we should point out that throughout 
the calculation we will omit the collisions. We will discuss the impact of collisions on the 
dispersion relation when we apply the results to the case of Farley-Byneman turbulence 
using the relaxation model for the collision operator already shown above. 

We now separate the Green's function into two parts; an unperturbed part and a 
perturbed one i.e., 

g = g(°) + gW (2.16) 

satisfying the following equations 

^ + iC (x,v,t)) a( )(x,v,t;x , ,v / ,O = ^(x-x , )<J(v-v / ^(t-O- (2-17) 

The solution to equation ( |2.17| ) can be written in the following formal form 

g<® (x, v, t; x', v', = 9(t - f )*(x(t) - x'(t'))£(v(t) - v'(tO), (2.18) 



where Q(t — t') is the Heavyside step function. Substracting equation (2.17) from equa- 
tion ( |2.14| ) leads to the following expression for the perturbed part of the Green's function 

e (1) (x,v,t;x',v',0 

= -iJ*<&'J d*" J v , ^(°)(x,v,t;x /, ,v ,/ ,t ,/ )A(x",v",0^(x ,/ ,v ,/ ,t";x , ,v , ,t / ). 
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Define the average Green's function G(x, v, t;t') as follows 

G(x,v,t;0 = J dx' J dv' (g(x,v,t;x'y,t')), 

where (• • •) represent an ensemble average. Taking the ensemble average of equation ( p. 14 ) 
and integrating over x' and v' leads to the following equation 



dG(x,v,t;t') 



+ J dx J dv' (iC(x,v,t)g(x,v,t;x' ,v',t')) = 5(t -t'). (2.19) 



dt 

Substituting equations ( ^T5|) , (|2~T^ ), ([D||) and ( |2~T7|) into equation ( pj|) we obtain 



dG ^ t ') + J & [ dV t dt" I dx" I dv" 

x (A(x, v, t)Q^ (x, v, t; x'W')A (*", v", t")0(*", v", t"; x', v', t'))- 



(2.20) 



Following the work of Ishihara et al. [12] and references therein we can simplify equa- 
tion ( |2.20| ) after making the following approximations in order to evaluate the second 
term on the left hand side. We assume that the process under consideration is Gaussian, 
which leads to 

/ dx' [ dv' (£ 1 (x,v,t)e (0) (x,v,t;x",v",t")A(x /, ,v ,/ ,t")e(x",v",t ,/ ;x / ,v / ,t')) 
J J (2.21) 

= (A (x, v, t)g^ (x, v, t- x",v", OA(x", v", t"))G(x", v", t"; t'). 

The direct interaction approximation makes a further approximation and replaces in 
the right hand side of equation ( |2.21| ) by Q. Finally, the last approximation consists of 
replacing G(x", v",t";t') by G(x, v,t;t'), and £i(x", v", t") by L x (x,v,t") to obtain 

dt Jt> J J (2.22) 

x (d(x, v, t)G(x, v, t; x" , v", i")>Ci(x, v, t"))G(x, v, i; t') = 0. 

Noting that C\(x, v,t) is given by equation (|2.15| ) we can rewrite equation ( |2.22|) as a 
diffusion equation. However before we do this we will assume that the ensemble aver- 
aged distribution function changes very slowly through secular changes of the integrals 
(constants) of the unperturbed motion, i.e., we assume that 

(/e(x,V,t)) = (f e (x,V±,Vpt)), 

where x represent the guiding center coordinates of the electrons and v±, vn the per- 
pendicular and parallel velocity components, respectively. This assumption allows us to 
express the velocity gradient in terms of gradients in the guiding center coordinates. The 
xpression for iC\ is as follows 

<5E k (t) • V v = -i$ k (t)k ■ V v 

which requires 

k . Vv = k ± ■ v±— /- + *„ /- + ^ • V, (2.23) 

V±OV± " OVil il e 
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where V represents the gradient with respect to the guiding center coordinates x. Note 
that the last term on the right hand side of equation ( 2.23| ) comes from the x B 
contribution. The diffusion equation can then be written in a compact form as follows 



dG(X,t;t') _ ^ d d 

at " \ dx~ a af5 ~dx~^ 



where X represents the guiding center coordinates as well as the perpendicular and parallel 
velocity components, i.e., 

X = {x,^,^} 

and where by definition the diffusion coefficients are given by 

D af) = ~(AX a AXp) 
more explicitly the diffusion coefficients are given by (when setting t' = to = 0) 

f>±± = ~(85i x (t)6x: ± (t)) 

= f Q dt" J dx" J dV (v B (x ) t)g(x ) v,tix",v",Ovs(x,0) ) 

= __!_ I dt" [ dx" I dV (v s (x,t)e?(x,v,t;x",v ,/ ,t ,, )<5E||(x,t ,/ )}, 

m e Jo J J 

= ^2 f dt " I dx " / dv " ^E||(x,t)a(x,v,t;x",v",t / ')5E||(x,t")}. 
m e JO J J 

The other diffusion coefficients, D„ x „ x , T) v±v .., T) v± ± can be expressed in a similar way. 
One can write a formal solution to the Green's function equation as follows 

g(x, v, t; x', v', t') = @{t - t')S(x - x')5(v - v')U(t, t'), 

where the time propagator U(t, t') is given by 

U{t, t') = exp f-i J £(x, v, t") dt"^j (2.24) 
This allows us to simplify the expressions for the diffusion coefficients 
I>±± = f dt" {v E (x,t)U(t,f)v E (x,lf% 

J 

D ± „ = - — [ dt" (v E (ic,t)U(t,f)SE^,lf% 

m e Jo 
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e 2 r* 

B V]n =— / dt"(<JE||(x,t)C7(t,t")*E||(x,t")>- 

m e JO 

The solution to the equation (|2.13| ) can be formally written in terms of the time prop- 
agator defined through equation ( [2.24| ), that is 

Sf e (x,v,t) = U(t,0)6f e (x,v,0)-i [ dt' ' tf(t,OA(/ e (x,«±,V|i,0>. 

Jo 

This leads to the following Fourier component of the fluctuating part of the distribution 
function 

5Wv) = -iCitoifefrv^vpt)) f dt'e^'H e- ik - x J7(i,t')e ik - x ), (2.25) 

where we have pulled out the velocity derivative of the averaged distribution assuming 
that the time dependence of the average distribution function is much slower than the 
time dependence of the orbits, and neglected the initial condition. 

Notice that the solution presented through equation ( 2.25| ) has imbeded in it two time 
scales. A fast time scale (associated with u>, i.e., fast oscillation time scale), and a slow 
time scale that should be related to the growth time scale. The slow time scale is the 
classical quasi-linear time scale required for saturation. The absence of a slow time scales 
leads to the absence of saturation. The Fourier transform is clearly a transform over the 
fast time scale. 

At this point we need to use a fundamental property of the time propagator, which is 
(U(t,t')g(x,v,t)) = g[x'(t')y(t')) 

and 

U(t,t') = exp dt" C(t")\ 

= exp dt" C(t")j exp ^-i dt" C(f)J = L/~ 1 (* / ) i )l7(Mo)- 

Substituting this result into equation (|2.25| ) leads to the following result 

(e- ik ' x [/(My k - x ) = (e- ik - x C/- 1 (*',*o)^(i,to)e ik - x ) = ( e - l ' k - x (*') e *-x(t)^ 

where x(i) and x(t') represent the exact orbits of the electrons. These orbits can be 
decomposed into the unperturbed orbits Xo(t) plus the perturbation due to the nonlinear 
effects of the random "bath" of waves Sx(t), i.e., 

x(t) =xo(t) +<5x(t), 

where the perturbation can be expressed as follows 

5v± 



Sx(t) = Sx(t) + <5xi| + 
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This finally leads to the following expression for the perturbed distribution function 

<5/eku,(v) = -i — Qto (k ± ■ V ± — + fc„/- + ^ • V) (/ e ) 

m e V Vj_ cwj. 11 OV\\ il e J 



(2.26) 

x / dt e (Mt-t')-*-(xo(t)-xo(t')))( exp ^ k . Ax(M'))>, 



where 

Ax(M') = <5x(t) - <5x(i')- 

At this point one can use the cumulant expansion (see for example Weinstock [20]) to 
write 



(exp (ik • Ax(t,f))> = exp ( f; l([ik ■ Ax(t, f)] n ) c j , 

\n=l ' / 



where (• • -) c is the cummulant, and note that this expansion can be reduced to a single 
term in the case of a random variable of a Gaussian process, that is 

(exp (ik • Ax(M'))) = exp (-±<[k ■ Ax(M')] 2 )) • (2.27) 

The exponent can now be expanded in the following form 

<[k • Ax(M')] 2 ) = k\ {{Sx\{t)) + (fej(0) - 2(8x\\{t)8x\\ {£))) 

+k ± • {(6Zl(t)) + {5St 2 ± (t')) - 2<$x L (t)$x L (t / )» ■ k ± 
+2k\\ ((6x\\(t)8x. ± (t)) + (Sx^SSc^t')) - {6x\\(t)5x ± tf)) 
-{Sx\\(lf)5Si ± (t))) -k ± 



+k --{-^ + ^ 2 a* )^ 

, 9 , f (Si ± (t)6v ± (t)) (SSt ± (t')5v ± (t')) {5St ± (t)5v ± (t')) 
+2k ^ • I ^ + ^ ^ 

n e )^ 

+ ^ / <gy ± (t)fa||(t)) + <«Vj.(Ofa||(f)> (*Vj.(t)fe||(0) 



(2.28) 



(Sv ± (t')Sx ll (t)) 



e 



fen. 



where the first two terms and the fourth term on the right hand side represent the par- 
allel and perpendicular correlation functions, while the other terms represent the cross 
correlation terms. 
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The diffusion coefficients are expressed in terms of the correlation function of the 
corresponding random forces due to the background oscillations. They can also be defined 
more accurately as the rate of time change of velocity variance around the mean value. 
The components of the velocity diffusion tensor can be written in the form: 

D VaV P = \j t {fo a {t)5Vp{t)), 

where the subscripts a and f3 represent the parallel and perpendicular components. The 
mean value of the parallel velocity is zero, while the mean value of the perpendicular 
velocity is given by the E x B drift due to the electric field of the electrojet. 
The random velocity components on the other hand are given by 



(t) = -— /"* dt' 6E\\(x(1/),1/), 



'o 

<Svj_(t) = --|cJE(x(t),t)xz, 

where we have set the initial time to be to = 0. We then have an expression for the guiding 
center spatial displacements 



Sx»(t) = -— f dh [ * dt 2 6E«(x(t 2 ),t 2 ), 

m e JO Jo 
C /"* 

6*±(t) = -- J dh (5E(x(ti),ti) x z). 



Parallel Velocity Diffusion. Following the work of Ishihara et al. [12], the parallel 
diffusion coefficient D v ^ v ^ can be obtained by evaluating the parallel velocity correlation 
function 

(5vl(t)) = 4 f dt' f dt" (<5£||(x(0,^£||(x(O,i")> 

m e JO JO 
2 rt r s 

= ^2 ds ds' {8E ll (yi(s),s)dE ll (x(s-s'),s-s')}. 

m e JO Js-t 

Using the definition of the exact orbits leads to 

o2 rt 



JO Js-t ,, 



ml 



It is clear from this equation that one needs the correlation of the spatial displacements 
to evaluate the parallel as well as the perpendicular and cross diffusion coefficients. 
In general, we need three correlations functions: {5x\\ (t)5x\\ (t')), (5x\\ (t)<5xj_(t / )), and 
(5xj_(t)(5x_|_(t / )). The first correlation function is given by 

{5x\\{t)5x\\{t')) = ^ f T dh [ ll dt 2 f T dt[ f l dt' 2 (5E\\(x(t 2 ),t 2 )SE\\(x(t> 2 ),t' 2 )) 

m e JO JO JO JO 
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following the arguments of Ishihara et al. [12], one obtains the following result for the 
parallel correlation function 



{Sx\\(t)Sx\i(t')) 



( -D m t' 2 (3t-t') for r>r', 
3 

-D m t 2 (3t'-t) for t'>t. 
. 3 



Perpendicular Guiding Center Spatial Diffusion. We now proceed to evaluate the 
perpendicular correlation function 

(5x x (t)Sx ± (t')) = F dt i r dt * ((<JE(x(ti),ti) x z)JE(x(i 2 ,i 2 ) x z)}. 
B Jo Jo 

It is important to note that this correlation function was approximated by Dum and 
Dupree [7] by 

((5x ± (r)5x ± (r / )}Dum-Du P ree OC D±±(t - r') 

using a Markovian approximation, and assuming the diffusion coefficient to be indepen- 
dent of time. However, such an approximation has its physical limitations as shown by 
Salat [18] and Ishihara et al. [12]. 

We now proceed to evaluate the perpendicular correlation function by noting that by 
definition the diffusion coefficient describing the diffusion of the electron guiding centers 
is given 

DjjL(t) = ~<« ± (t)«j_(f)) 

= \^2j t J o dt lj Q *2 <(£E(x(ti),ti) xz)(<5E(x(i 2 ),i 2 ) xz)) (2.29) 



c 2 







B 2 

This in turn leads, in the case r > r', to 

„2 



dh ((<5E(x(ii),£i) x z)5B(x(t),t) x z)). 



(<5x±(r)5x J .(r / )) = jf dti jf di 2 <(*E(x(ti),ti) x z)(<5E(x(i 2 ),t 2 ) x z)) 
c 2 f T r' 

+ 52 y, dtl _/ ^ 2 <(<5E(x(ti),i!) x z)(<JE(x(t 2 ),t 2 ) x z)). 
The first term can be approximated using the definition of the diffusion coefficient 
(5x ± (r)<5x ± (r / )} = 2D ±± (r')r' 

„2 



C 
+ ^2 



/ rfti / dt 2 ((<JE(x(ti),ti) x z)(5E(x(t 2 ),t 2 ) x z)). 

ir' JO 



Using the definition of the diffusion coefficient given by equation ( 2.2S| ) and changing 
variables allows us to rewrite the integral in the following form 

„2 



(£x ± (t)<5x ± (t')) = 2B ±± (t')t' + iL ' d£ dr} 



(2.30) 

x <(<5E(x(£ + r'), e + t 7 ) x z)(5E(x(r?), 7?) x z)). 
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At this stage we would like to consider the case where r r', i.e., we shall consider 
the case where \t — t'\/t' is a small parameter. This assumption is critical to the rest of 
the calculations in this paper. It is also important to note that in general r and r' are 
independent variables, which makes our assumption a critical one indeed. We then can 
write 

5E (x(£ + t'U + r') = <5E(x(r'), t') + ^JE(x(t'), r') ; 
where the convective derivative is defined by 



A 

(It 



7-(r7+vM.v). 



The first term when substituted into equation ( |2.30| ) leads to Dj_j_ (r')(r— r') and therefore 
the expression becomes 

(55L x (t)6Z ± (t j )) =D±±(r')(r + r') 

2 .r-r' r' ( d \ ( 2 ' 31 ) 

+^2 y o y o ^ <( ^e^t'),/) x z j (SE^r?) x z )>. 

We now use the following property of the diffusion coefficient 



r D ±± (t) = -j<(<JE(x(t),t) x z) (<fE(x(t),t) x z)> 



(2.32) 

C / " . .Id . . . \ 

+- 



B 2 



J dh (f^tfE(x(t),t) x (<SE(xfa) } fi) x z)). 



Consequently we can express the two-time correlation function by substituting expres- 
sion ( |2.32|) into expression ( |2.31| ) 

(5x ± (r)<5x ± (r / )) = D ±± (r')(r + r') + (r ~ ^ (^D^r') 
- ^((5E(x(r'),r') x z)(5E(x(r'),r') x z)> 

Therefore 

((5x ± (r)5xj.(r)) + (<5x ± (t')5x ± (t')> - 2(<5x ± (t)5x ± (t')> 
= 2D ±± (r)r + 2D ±± (r')r' - 2D ±± (r')(r + r') 

-<T ~ T 'f (^7 6 ^( T ') " ^((^(x(r'),r') x z)(<5E(x(t'),t') x z))) . 

The first four terms on the right hand side of expression ( |2.33j ) can be rewritten, after 
expanding Dj_^(r) around r', in the following form, 

2 (Dll(t) - D ±± (r')) r - (r - rf^Dn^) = (r 2 - r'2) ^Dii^), (2.34) 



(2.33) 



452 



A.M. Hamza 



where we have used 

D ±± (r) = T)±±(r') + (r - T ') -^jB ±± (r') . 
Finally equation ( |2.33| ) can be written as 

- S/Mr , MX* , (t')\ = (t 2 - t' 2 \ 

dr 



<<yxi(r)«i(r)) + (<5x ± (r')5x ± (r')) - 2(Jx ± (t),5x ± (t / )> = (r 2 - r'^D^^) 



,2.35) 

+ (r-r') 2 | I ((<5E(x(r'),r') x z)(ffi(x(r'),r') X z)>. 

The second term can easily be identified with (5vj_(t')<5v_l(t')), and can therefore be 
identified with a velocity space diffusion coefficient, which decribes random changes of 
the gyroradius and phase angle, and can therefore be written as 

(<5v ± (t')5v ± (t')) = 2D V±V± (t')t'. (2.36) 

Substituting expression (|2.36| ) into equation ( |2.35[) leads to 

(<5x l (t)5x_ l (t)) + (<5x x (t')<5x ± (t')) - 2(5x x (t)<5x ± (t / )) 

rfV2 /2x ^ / /n on 2/ ^/A2^/ Z' A^u±(0 A (2.37) 
= ( T ~ T )^7 D ^M T ) + 2SJ e( T T ) T I ^2 ] • 

One can finally use the approximation of Dum and Dupree [7], that is the velocity diffusion 
makes approximately an equal contribution to the diffusion of the guiding centers, i.e., 

D v±vAr')\ ~d ±±(t ') (2.38) 



which finally leads to the expression for the perpendicular correlation functions ( |2.37| ) 

(55c ± (t)6Z ± (t)) + (S5t ± (r')SSt ± (T')) - 2<<5x ± (t)5x ± (t')) 

r (2.39) 
= (t 2 - t' 2 )^T) ±± (t') + 2tf e (r - t') 2 t'D ±± (t'). 

In the case r' > r we just exchange the role of r and r' in the expression ( |2.39| ). 
Perpendicular Velocity Diffusion. For the velocity diffusion, with r > r' we need to 
evaluate the following expression 

<<5v ± (r)<rv ± (r')) * (^v ± (r')^(r')) + ^^^7(*v ± (t')<5v ± (t')>. (2.40) 

Using the definition of the velocity diffusion we can express the right hand side of equa- 
tion ( p. 40 ) as follows 

(6v ± (r)6v ± (r')) ~ 2B v±v± (t>)t> + (r' - T) ^^)ll 

which finally leads to the contribution of perpendicular velocity diffusion 
(6v ± (t)6v ± (t)) + {8v ± (r')Sv ± (r')) - 2(8v ± (t)5v ± (t')) 

~ 2B v±v± (r)r - 2B v±v± (t')t' - 2(r - r')B v±v± (r') - 2r'(r - T >)^J^Jzl . 
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Expanding the first term on the right hand side leads to 

D v±v± (r) - D^(r') + (r - r' f^j/^ 
which leads to 

(8v ± (t)5v ± (t)) (5v ± (t>)5v ± (t>)} {SvJ^dvJj^ ^ _ f , 2 dt) ±± (r') 

+ w e ~ Z{T T) dr' ■ 



where we have used equation ( 2.38 ). 

Parallel Velocity- Perpendicular Spatial Cross Diffusion. We have now evaluated 
the parallel and perpendicular correlation functions, we therefore only need to calculate 
the cross correlation, which is given by 

<<yx.|(T)<5xjV)) = —4 [ T dh [ h dt 2 f dt[ (8E ll (x(t 2 ),t 2 )5B(x(t' 1 ),t> 1 )xz)).(2Al) 

We proceed in a similar fashion to the previous cases, and start by defining the cross 
diffusion coefficient 



1 e c 

1 e c 
2ml B 



t 

dti (6E\\(x(ibi),t 1 )5E(x.(t),t) x z)) 

t 

dti (<5S||(x(t),i)<5E(x(i),i) x z)). 



In the case r > r' we can decompose expression (2.41) and write it in the following form 

(fo|,(i)<5x ± (i)) = \ — ^ f dh f dt 2 {SEn{x{h),h)SE{x(t 2 ),t 2 ) x z)} 
l m e d Jo Jo 

which allows us to approximate the integrant of equation ( |2.41| ) as follows 
— t dt * f dt[ <5£7|,(x(t 2 ),t 2 )$E(x(f 1 ),f 1 ) xz)) 

pmin(ti ,t' ) />min(ti ,r') 

/ dti/ dti (^||(x(t 2 ),t 2 )«5E(x(t , 1 ),t' 1 ) xz)}. 



e c 
Trip 73 



The right hand side of this equation can easily be identified as 

(£*U||(£)<5x ± (t)} t=min(tl)T ,) « 2D V]i± mm(t 1 ,T f ). 
Expression ( |2.41[ ) can now be approximated as 

(Wr^XiCr')) = f dh 2D ± h + f dh 2D , r' = 7J ± (2r - tV- 
On the other hand 

(Wt')£x ± (t)) = — -| T rfti Tdt, [ T dt[ (<5£,|(x(t 2 ),t 2 ) (JECx^),^) x z)>. 

m e -D Jo JO JO 
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This expression can be evaluated for r > r' following the same procedure as before to 
obtain 



{Sxu(t')Sx±(t)} = [ dti 2D V , min(t 

Jo 11 



,r) = A,.. ,r' 2 . 



This finally allows to sum all the cross correlation functions that appear in the expres- 
sion ( §jH) 



(<yxii(t)5xj.(t)) + (<yx.|(0<5x ± (0> 



<<$sii(t)<$x ± (t / )) + (fey (0^±(*)> + ,±(0(* - + (* - *0< 



<9A 



— v 



at 1 



Cross Perpendicular Velocity- Space Diffusion. One can easily show that the cross 
perpendicular diffusion (velocity-space) can be written as follows 



2 dT> v±± (T>) 

8t' 



(5v ± (t)5x ± {t)) + {6v ± (t')SSL ± (t')) 

= (5v ± {T)5St ± (T')) + (Sv ± (T')5St ± (T)) + 2(r - t) 

Cross Velocity diffusion 

(6 Xll (r)8v ± (r , )}=D Vl[Vx (T , y(2r-T% 
(6x ll (r')6v ± (r))=D VllV± (r')r , \ 

which finally leads to the following expression 

(fe||(t)5v ± (t)) + (fe||(t')5v ± (t / )) 

3D 

= (fo||(t)*v ± (f )> + (S Xll (t')Sv ± (t)) + D Vl{Vx (t')(t - t'f + (t- t')t 2 —^. 

We now are able to express equation ( 2.34| ) in terms of components of the diffusion 
tensor, that is for t > t' 



k • Ax(t, t')] ) = kU -D vn (t 3 + t' 2 (2t' - 3i)) + k ± 



d 







+ 2(t-t') 2 —\t) ±± (t') + 



D u±± (i') 



+ 2^(t-t')1 / D il (t') k 



(2.42) 



(< - 



/\2 



£„„±(0 + 



k, . 
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At this stage we need to solve self-consistently for the diffusion coefficients introduced 
above. First, we start with the parallel diffusion 

1 d 



D 



2 dt m 2 

\_d_e^ rt rO 

2dt ml 



+ 



1_ d_^_ 

2dtml 



[ ds f ds / V|(5S[| k | 2 e i ^ k -t xo W- xo ('-* , W- Wk "')(e ik - Ax ("''-' / )> 

Jo Js-t k 

r t ,o (2-43) 
/ ds / A/Vl^HklV^ 1 

JO Js-t k 

f ds I ds'V|(5£|| k |V 
o Jo . 



•(x (s)-x (s— s'))— u; k s') ^ik-Ax(s,s-s')^ 



3 i(k-(x (s)-x (s-s'))-Lj k s / ) ^ik-Ax(s,s-s') 



Note that different expressions ought to be used for Ax(s,s — s') for each of the last 
two integrals of expression ( 2.43 ). Then using equation ( [2.27D and defining the following 
diffusion coefficient we obtain 



V ±1 _ = + 



a, ' 



„ — f II V I 

Case when s' > 



exp (-~([k • Ax(s,s- s')] 2 )) = exp | -yA, n s' 2 (3s - 2s') 



/2~ ( /\ | 2 ,^2||±( S -S' 

s S|i(s — s ) + s s - 



ds 



x exp < — k 



fi^ ,2 ( S _ s ')t) ±± (s - s') + s' ( 8 - ^ ) ^D ±± ( S - *') 



s'\ d 



(2.44) 



2 / ds 



i2 



dV ±± (s - s') 



Os 



•k 



Case when s' < 

exp 



1 \ f k 2 , 

-([k- Ax(s,s- s')] 2 )J = exp < -yA, ||t)|| s /2 (3s - s') 



s' 2 V ll± (s)-s'(s-s') 



x exp | — kx • 

,2 3£j_l(*) 



s'\ d 



(2.45) 



'sDn(s) - s' [ s - - J ^D ±± (a) 



(9 .5' 
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Using the last two expressions ( 2.44 ) and Q2.45 ), and changing variables in the second 
integral, we obtain the following expression for the parallel diffusion coefficient 

e 2 x f* 

Dv \\ v \\ = ~2 z2 l^llkl 2 / ds COS ( k ' ( X o(*) ~ x o(* - S)) - UJ k s) 
m e v. JO 



x exp < ——D v ^ v ^ s 2 (3t — 2s) — k\\ 



s 2 Vu ± (t-s) + t 2 s- 



dt 



x exp < — k± 



s s d ~ 



nis 2 (t - a)T> ±± (t - a) + s(t - -)irD±±(* - s) 



2 '(it 



2 dV ±± (t-s) 



at 



•k 



A similar calculation for the perpendicular components of the diffusion tensor as well as 
the cross components, leads to 



c 2 f* 
Djj.(t) = ^2 ^(^ E k x z)(^ E k x z) / cos (k • (xo(t) - xo(t - a)) - w k s) 

-D v ,, v ,,s 2 (3t - 2s) - fcii 



x exp 



k 

fey 

~3 



_ dVn ± (t - s) 
s 2 V\\ ± {t-s) + t 2 s-~ 



Ot 



x exp 



s ^ d ~ 



n 2 e s 2 (t - a)D ±± (t - a) + s(t - -)-D ±± (t 



2 SP ±± (i-a) 



and similar expressions for the cross diffusion coefficients. 

It is clear from these results that the resonance broadening effects due to scattering of 
electrons by the Modified Two Stream Farley-Buneman waves can not be accounted for by 
just replacing uj by uj + ik^D* in the resonant part of the dispersion relation. It is obvious 
from the expressions above that there is a complex time dependence of the diffusion 
coefficients. Moreover, most of the published work (Dum and Dupree [7], Sudan [19], 
Robinson [16], Robinson and Honary [17]) ignores the cross correlation and consequently 
the cross diffusion. 



3 The Nonlinear Dielectric Function 



In order to obtain the dielectric function one has to use the results of the previous section. 
The expression for the fluctuating part of the electron distribution function ( [2.26 ) can 
now be written in the following form after using equation ( 2.42j ) 
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£/eku;(v) = -*Aku;(/e) 

x jf dt' e (i»(t-t')-iMMt)-Mt>))) exp (f + t' 2 ( 2 £' - 3t)) I 

x exp |-k ± • ^L-^-D^f ) + (t - t'?^ 1 (3.1) 

+ n 2 e (t - t') 2 t'i) ±± (t')) ■ k ± - k\\ (v\ ]± (t')(t - t'f 
+ t\t-t' 



dt' 



using this expression (|3.1|) along with equation (|2.3|) in Poisson's equation, and changing 
the integration variable from t' to s = t — t 1 leads to the nonlinear dielectric function 



e(k,w) = 1 — — : I dv(k-v_L — — — + k\. 

If Li " 

+ ^ " V_l) (/ e > e (--*-(x (t)-xo(t- S ))) exp \_ k l DvnS 2 {2t _ 2s) 



X 



expj-^ [s%^t-s) + t 2 s d -^ t S) ^j .k ± | (3-2) 
x exp |-k± • fs(t - |)^D_l±(* - s) + ^(t _ a )D_L ± (t - s ) 

This dispersion relation differs considerably from that used by Robinson and Honary 
[17] and Robinson [16] in many aspects, and therefore the consequences on the physics of 
irregularities in the auroral as well as equatorial E regions are significantly different. We 
will show in the next section how the accurate resonance broadening calculation affects the 
threshold for Farley-Buneman instablity. We will also show how it affects the important 
problem of aspect angles. 

In order to extract the information hidden in the dispersion relation we have to evaluate 
the time integral in the expression (|3.2j), a nontrivial calculation. 

A final note on the time dependence. As mentioned earlier in the paper, the time 
dependence that appears in the right hand side of the dispersion relation is a slow time 
dependence necessary for energy and momentum balance. In other words, the time depen- 
dence is necessary for wave saturation. This is a classical problem. Linear theory ignores 
the slow time dependence and predicts a time independent growth rate which in turn sug- 
gests that waves will grow indefinetly. Quasi-linear theory remedies this critical problem 
by introducing a slow time dependence in the background distribution function. We have 
retained both the fast and the slow time dependence and Fourier transformed over the fast 
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time scale. The slow time scale is associate with the diffusion time scale in the classical 
quasi-linear theory. This same time scale reappears in the resonance broadening analysis. 



4 The Farley-Buneman Case 

In order to simplify the results and obtain a direct comparison with the classical Farley- 
Buneman results we will assume that SE^ = 0, and the we have isotropic turbulence. 
This allows us to eliminate the parallel and cross diffusion effects, and use a simplified 
dispersion relation 

uli wjL f , / Id k x z 



e(k, u) = 1 - j2T pt , z-f / dv ( k • v ± — + ^ • V ± ) (/ e 

- k 2 ±- + iv in u) k 2 J V v ± 9v± n e 



x jTds e (^ s -ik.(x (t)-x (t- s ))) exp |_ fe 2 ^ _ |)^Z>*(t - s) 
+ ^ s 2 (t - s)D*(t -s) + s>P^ ~ S A 1 



where we have assumed an isotropic spectrum for simplicity, and replaced Dj_j_ by D*. 
Then assuming a slow time dependence of the diffusion coefficient D* we can further 
simplify the expression for the dielectric function to obtain 

wji w 2 e f ( 1 3 kxz 

/ dv(k-v ± — — + — ^ 



e(k, cu) = 1 - , / civ k • v ± — — + -— • V ± (/ e 



x r ds e (i««-*-(xo(t)-xo(t-«))) {_jfe2 n 2 s 2^ _ s j_ D *^ t _ g )| 
JO 

At this point one can explicitly express the unperturbet orbits of the electrons in terms 
of Bessel functions, that is 



+oo 



e (ius-iMMt)-Mts))) = J- Jn(k±p±)J m (k±p±) 



TL,m= — CO 



x exp |i (to — k • — nVt^j s + i(n — m)(f) — i{n — m)fi e i j 

the integral over the velocity allows us to reduce the double sum to a single sum by 
integrating over c/> to obtain 

to 2 to 2 f +oc 



n=—oo 



where the resonance function R is given by 



R(u - k ■ - nVL e ) 

rt (4.1) 
= J ds exp^ifu-k-v^ -nn^js-kln 2 e s 2 (t-s)D*(t-s)y 
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It is clear from the expression ( |4.l| ) that the resonance function is similar to that of the 
unmagnetized case derived by Ishihara et al. [12] and Salat [18]. Therefore one can extract 
similar properties of the time integral and consequently obtain some useful information 
regarding the broadening of the wave-particle resonance. If we were to make the same 
change of variables 

t k = (k 2 ± n 2 e D*)-^, T=—, U = u - k • v ( F 0) - nfi e , R = — 

T~K T K 

then the resonance function can be written as 

rT 



R(U,T)=[ dS exp{iUS - S 2 (T - S)}. 
Jo 



10 

Note that only the real part of R appears in the expressions for the diffusion tensor 
components. 

The final question to be addressed regarding this problem is that related to collisions 
between electrons and neutrals. It has been established through linear theory, Kadomt- 
sev [13], Coppi and Rosenbluth [2] and Hendel et al. [9], that collisional effects enter the 
dispersion relation through the introduction of a collisional damping d c ° 11 depending on 
the wave frequency, and which can be expressed as follows for electron neutral collisions 

icoll _ u e _ 1,2 n co11 

2[y en - luj) 

where D c ° 11 is the collisional diffusion coefficient. 

When combining the last two equations, we obtain the correct dispersion relation for 
Farley-Buneman waves. To be more specific we can integrate the Vlasov equation with 
the Bhatnagar collision operator, one can formally solve for the distribution function and 
therefore the charge dewnsity to obtain the dispersion relation. The steps of this procedure 
are described below. 

We start with the collisional Vlasov equation, which we write in the following form 

^JT + "endfe = —V<$> ■ V v / e0 + V en —8n(x, t), 

at m e uq 

where the total time derivative d/dt is the derivative along the perturbed particle orbits. 
The solution to this equation can be written in an integral form 

Sfe = /eo(v)e-"™ ( 

+v en e-'*» t fdt' { e —V^(t'),t') ■ V v / e0 (v) + M^ n ( x (i') ;t ')\ e -«W 

Jo { m e v en n J 

writing 

5f e (x,v,t) = 5f e]iLU exp{z [k • x(i) -wt]}, $(x,t) = exp{i [k • x(t) - tot]} 
leads to 

Sfe-kw = feO exp{-i [k • x(t) ~(U + lVen)t}} 



+ 



f dt (i^ / e0 k • V + A 5nk \ ^«^n)(«') {e -*[k-W*)-^))] ), ' 

Jo [ T e n J 
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where x(t) = xo + 5x(t), with xo representing the unperturbed orbits, while 5x is the 
perturbation due to the random electric fields. Note that we have assumed that f e o is a 
drifting Maxwellian. Neglecting the first term in equation fl4.2| ) for long times t —* oo, 
and using the results of the previous section we obtain 

Sfeku = [ ds li^^feok ■ v + v en —5n VuJ 

JO { T e n 

xe i((u J +iv en )(t-t>)-MMt)-x Q (t-s))) exp f-i([k ■ Ax(i,i - s)] 2 ) 

Using the results of the previous section where the Resonance function was calculated, we 
can deduce the expression for the fluctuating part of the distribution function, and then 
integrate over velocity to obtain an expression for the density 



5n ku; ,e$ kt(J f J /"*_,/, 
= i „ 



d v / ds — k- ve i((^+^e»)(*-*')-k-(xo(t)-xo(t-s))) 



™0 J-e J Jo n 

xexpf-i([k- Ax(M-s)f) 



dx ^ d S Il e i{(u-Hven)(t-t')-k.(Mt)-Mt-s))) 



+v e J-^ L [ dv [ ds — 
no J Jo n 

xexp f-^([k- Ax(M-s)f) 

Substituting this expression along with the expression for the ion charge density in Pois- 
son's equation we obtain the dispersion relation 



io 2 -k 2 ^-+ iv in io 

J JO nn 



^ 2 A|) e ^ _ f (7 V f* ds A e <(^+^)(*-*')-k-(xo(t)-xo(t-s))-|<[k-Ax(t,t- S )] 2 )) N 
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5 Summary 

We have shown that the components of the electron diffusion coefficient are time depen- 
dent and the conventional result that suggested replacing a; by a; + ik\D* in the resonant 
part of dielectric function is not valid. This alters a number of results obtained through 
the application of the classical resonance broadening calculations of Dum and Dupree [7] , 
such as Sudan's results [19] and Robinson's results [17] concerning the thresholds of the 
Farley-Buneman and Gradient drift instabilities in the ionosphere. The correct results 
will be presented in a subsequent paper to be submitted in the near future. We have 
aslo added the parallel diffusion as well as the cross diffusion coefficients. The former is 
identical to the one dimensional analog derived by Ishihara et al. [12] and Salat [18], the 
latter has never been calculated explicitly. 



Resonance Broadening Theory of Farley-Buneman Turbulence 



461 



Finally, we have explicitely derived the dispersion relation for the Farley-Buneman 
waves using the improved resonance broadening formalism. Further details on the Farley- 
Buneman thresholds and transport will be published in the near future. 
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